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I .  INTRODUCTION 


Recent  efforts  here  in  the  field  of  Quantum  Optics  have  led  to  the 
generation  and  establishment  of  important  results  and  findings  in  the 
areas  of  swept-gain  superradiance  [1,2], laser  induced  chemistry  [3] 
and  optical  bistability  [A],  These  areas  constitute  examples  of 
collective  phenomena  in  light-matter  interactions  which  have  incurred 
widespread  interest  and  theoretical  and  experimental  activity,  not 
only  from  the  standpoint  of  the  fundamental  understanding  of 
cooperative  effects,  but  also  from  the  motivation  stemming  from 
obvious  potential  applications. 

State-of-the-art  results  in  the  areas  of  superradiance  and 
cooperative  effects  were  first  reported  here  in  December  1976  as  part 
of  the  US  Army  sponsored  Symposium  on  New  Laser  Concepts. 

Subsequently,  these  presentations  were  published  by  Plenum  Press  as 
part  of  the  text  Cooperative  Effects  in  Matter  and  Radiation,  edited 
by  Bowden,  Howgate,  and  Robl  [5],  In  April  1977,  the  first  swept-gain 
superradiance  experiment  [1]  was  performed  at  our  laboratory  using 
C02~pumped  CH3F.  As  part  of  a  coordinated  effort,  our  theoretical 
program  has  modeled  superradiance  in  a  collection  of  three-level 
molecules,  where  the  dynamical  transient  effects  of  coherent  optical 
pumping  are  explicitly  treated  [6],  Here,  it  was  shown  that 
two-photon,  coherent,  stimulated  Raman  transitions  play  an  important 
role  in  the  evolution  of  macroscopic  collective  polarization.  The 
connection  between  pure  superradiance  and  swept-gain  superradiance  was 
then  explicitly  established,  both  experimentally  and  theoretically. 

Also,  the  influence  of  chemical  reactivity  and/or  dissociation  on 
laser-induced  population  transfer  in  multilevel  systems  has  been 
studied  [7]  and  simple  analytical  models  [3, 7, 8, 9]  have  been  devised 
to  determine  the  nature  of  essential  fundamental  processes  in  larger 
systems.  Conditions  under  which  population  can  be  trapped  in  excited 
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states  in  multilevel  systems  undergoing  dissociation  and/or  reaction 
have  been  established  [3].  Also,  a  convenient  and  economical  computer 
algorithm  has  been  developed  from  which  a  computer  simulation  model 
was  established;  and  rules  for  the  reduction  of  the  effective  manifold 
of  states  in  a  multilevel  laser  pumped  system  have  been  devised  which 
can  significantly  increase  the  efficiency  of  and  decrease  the  cost  of 
machine  computation  [9]. 

Likewise,  investigation  of  the  instabilities  in  the  ground  state 
of  the  Dicke  Hamiltonian  [10,11]  has  resulted  in  the  establishment  of 
optical  bistability  as  a  first-order  phase  transition  far  from 
thermodynamic  equilibrium  [4]  —  which  is  caused  by  virtual  photon 
exchange  in  atomic  pair  correlations  in  a  mean-field  interaction. 

This  observation  has  led  to  a  demonstration  of  the  precise  formal 
connection  between  optical  bistability,  superfluorescence,  super¬ 
radiance,  and  the  BCS  theory  of  superconductivity.  To  emphasize  our 
interest  and  support  in  this  dynamic  new  area,  an  International 
Conference  on  Optical  Bistability  is  being  sponsored  jointly  by  the  US 
Army  Research  Office  and  this  installation,  to  be  held  3-5  June  1980 
in  Asheville,  North  Carolina,  USA. 

It  is  not  the  purpose  of  this  paper  to  review  in  detail  our  total 
program  which  is  (on  the  whole)  well  documented  in  the  open 
literature,  but  rather  to  focus  attention  on  several  critical  areas  of 
thrust  and  current  momentum.  With  this  in  mind,  we  shall  attempt  in 
the  remainder  of  this  discussion  to  outline  and  present  some  of  our 
results  in  several  areas  of  current  endeavor.  Specifically,  Section 
II  will  consider  a  rigorous  procedure  for  systematically  handling  the 
general  interaction  of  a  molecular  gas  with  the  electromagnetic  field, 
while  Sections  III  and  IV  will  deal  more  directly  with  specific 
models  that  enhance  our  understanding  of  and  effective  treatment  of 
cooperative  phenomena  in  light-matter  interactions.  Section  V  will 


provide  a  brief  summary  of  our  conclusions  and  future  goals  and 
direction. 

II.  INTERACTION  OF  A  MOLECULAR  GAS  WITH  THE  ELECTROMAGNETIC 

FIELD 

In  recent  years  a  number  of  papers  [12]  have  considered 
theoretical  models  for  the  interaction  of  molecular  systems  with  the 
quantum  electromagnetic  field.  Problem  areas  centering  around  the 
treatment  of  the  molecular  system  Itself  (i.e.,  motion  of  the  center 
of  mass)  and  restrictions  pertaining  to  field  origin  and  structure, 
etc.  tend  to  impact  the  formalism  such  that  one  neither  knows 
implicitly  the  approximations  which  have  been  made  in  the  model  nor 
how  one  could  improve  on  the  results  of  the  calculations.  Ideally, 
the  basic  equations  of  the  theory  should  be  physically  “exact”  in  the 
sense  that  they  are  derived  from  first  principles;  they  should  be  of  a 
form  which  makes  the  physical  and  mathematical  significance  of  various 
approximations  reasonably  clear;  and  they  should  permit  any  desired 
physical  effect  to  be  calculated  with  sufficient  effort. 

In  order  to  achieve  some  of  these  goals  for  the  case  where  the 
molecular  system  is  a  gas,  a  complete  second-quantized  Hamiltonian  for 
a  system  of  nonrelativistic  massive  charged  particles  interacting 
electromagnetically  with  each  other  and  with  the  quantum 
electromagnetic  field  has  been  formulated  (by  C.A.C.)  subject  only  to 
three  a  priori  approximations  or  constraints:  The  nuclei  of  the 
molecules  can  be  considered  as  elementary  particles,  all  massive 
particles  can  be  treated  nonrelativistically,  and  wall  effects  on  the 
internal  states  of  the  molecules  can  be  neglected.  A  key  step  in  the 
formalism  of  the  theory  is  the  choice  of  an  appropriate  set  of 
operators  for  the  description  of  the  molecular  system.  For  this 
purpose  we  introduce  the  operators  P^CX, t)  which  destroy  a  molecule 
in  internal  state  k  with  renter-of-mass  position  X  at  time  t,  along 
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with  its  adjoint  p£(X,t)  which  creates  such  a  molecule.  A 
Hamiltonian  is  then  used  to  derive  equations  of  motion  for  the 
Pk(X,t)  which  contain  only  these  operators,  their  adjoints,  and  the 
electromagnetic  vector  potential.  The  current  term  in  the  equation 
for  the  electromagnetic  vector  potential  is  also  expressed  in  terms  of 
these  molecular  operators.  The  resulting  set  of  operator  equations  is 
closed,  forms  a  complete  basis  for  the  description  of  the  interaction, 
and  has  the  convenient  feature  that  in  it  the  internal  molecular 
states  (which  are  involved  in  electromagnetic  transitions)  and  the 
molecular  center-of-mass  coordinates  (which  play  a  role  in  transport, 
density  fluctuations,  Doppler  shifts,  etc.)  enter  in  different  ways. 
This  fact  is  advantageous  in  allowing  use  of  experimentally-determined 
parameters  in  numerical  calculations  for  specific  systems,  and  it 
permits  one  to  see  in  what  ways  various  approximations  in  treating  the 
molecular  system  interact  with  one  another. 

The  total  Hamiltonian  for  the  system  in  rationalized  Gaussian 
units,  nonrelativistic  in  the  massive  particles,  is 


where  ^m(2£»t)  is  the  annihilation  operator  for  a  particle  of  type  m 
(i.e.,  electron,  nuclei)  at  position  x  and  time  t,  A(x,t)  is  the 
electromagnetic  vector  potential  in  the  Coulomb  gauge,  qm  and  ^ 
is  the  charge  and  mass  of  the  mt'1  species  of  particles, 
respectively,  etc. 
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We  now  define 


Pk(I»t)  *  j$k(y)  n  d3y^ 

I  j 


relevant  to  a  molecule  where  £1  is  the  coordinate  of  its  ith 
particle,  X  is  the  coordinate  of  its  center  of  mass  (x.  *  £  +  £ 
b^^y^)  ,  4>k  denotes  the  k  eigenstate,  and  tp^  is  the  J 

annihilation  operator  of  its  i**1  particle  (to  be  distinguished  from 
i|/m,  as  m  denotes  specific  type  while  i  specifies  a  given  particle). 

It  is  obvious  from  this  definition  that  Pk(X»0  destroys  a  molecule 
with  internal  coordinate  k  and  center-of-mass  coordinates  X  at  time  t, 
and  correspondingly  that  P^(X, t)  creates  a  molecule  of  this  kind.  The 
field  operators  Pfc^k  are  convenient  in  studying  molecular 
interactions  with  the  electromagnetic  field  because  they  reflect 
properties  of  an  entire  molecule  rather  than  of  the  individual 
particles  that  constitute  the  molecule. 

We  now  develop  equations  of  motion  in  which  only  the  operators 
pk*pk  and  A  appear,  together  with  constants  and  functions 
determined  by  the  Internal  molecular  states.  This  set  of  three 
equations  now  forms  a  complete  and  closed  set  of  operator  equations 
for  the  description  of  the  resulting  N-molecule  system  interacting 
with  the  quantum  electromagnetic  field  —  the  only  approximations 
being  those  three  a  priori  constraints  specified  earlier.  For 
simplicity  an  additional  dipole  approximation  can  now  be  introduced  on 
the  assumption  that  for  states  of  interest  the  electromagnetic  vector 
potential  varies  slowly  enough  over  molecular  dimensions  that  this 
variation  can  be  neglected.  In  situations  requiring  greater  accuracy, 
this  (optional)  approximation  could  just  as  well  be  omitted  in  the 
resulting  calculation  of  the  dynamical  properties  of  the  system. 
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Our  basic  procedure  is  now  to  use  these  operator  equations  to 
generate  equations  for  a  hierarchy  of  Green  functions.  By  making  var¬ 
ious  factorization  assumptions  one  may  then  retain  as  much  of  the  ex¬ 
act  content  of  these  equations  as  is  necessary  and/or  feasible  in  the 
case  of  interest.  We  therefore  define  the  Green  functions 


Pkt(X,X')  *  b,  <  p£(x)  PA(X')  >  , 


a(X)  *  <  A(X)  >  , 


where  t  is  understood  and  bj  and  b2  are  normalization  constants. 
Under  a  rational  factorization  procedure. 


A 

V 

ii 

A 

V 

A 

l> 

V 

<  P^A2  >  = 

<  pfo  > 

A 

l*» 

V 

(4) 


The  resulting  Green  function  hierarchy  can  be  truncated  to  Include 
only  those  functions  defined  under  expression  (3).  In  this  minimal 
''semi-classical"  approximation  (that  includes  also  the  dipole  ap¬ 
proximation),  equations  of  motion  for  Pk£(X»X' » t)  and£(j(,t)  assume 
a  representation, 

•  [m  -  i'2>  -  <VEi>]  -  V-  2 

k,k,k0 
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-  ^I^jtr-i*)  \k1ic2it3<Ml'.*-.x,>] 

'  £  ^([“k'k  4'k.'i(-)Pk't(-,-)  *  “ik‘  4k' 

x  id-)  ■>„„. <*,*•)]  -fcSjrjp^x.x') 

x  [£2(X)  -  *2(X')]  . 

£(x.t)  -  £  i  [zv4tpk  t 


( 


XX  -  V2 
?  3? 


-  ^^(x)EPkk(x,x)]f 


(5) 


where  [ - ]t  is  the  solenoid  part  (zero  divergence)  of  the  vector 

quantity  in  brackets.  The  role  of  parameters  such  as  transition 
energy  w^k.' ,  associated  dipole  moment  d^k'x  intermolecular 
coupling  Qk.jk2k3k4»  etc*  is  obvious  upon  inspection  of  the 
Equations  (5),  which  do  not  form  a  closed  system  because  of  the 
presence  of  the  functions  Pkik2k3k4*  if  the  intermolecular 
interactions  are  to  be  treated,  one  must  either  include  additional 
Green  function  equations  or  make  additional  approximations  to  express 
the  Pk!k2k3k4  in  terms  of  the  P^jk^* 

Without  the  presence  of  the  Intermolecular  term,  Equations  (5)  do 
become  a  closed  set  describing  the  evolution  of  the  system.  They 
describe  the  processes  of  absorption  and  emission  of  radiation, 
Doppler  shifts,  density  fluctuations,  etc.,  but  contain  no  direct 
field-molecule  correlations  and  cannot  produce  the  spontaneous 
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emission  of  radiation.  To  proceed  to  the  next  higher  order  of 
accuracy  one  can  drop  the  assumption  that  <P[CP^A)  factorizes  and 
define  a  new  Green  function  proportional  to  this  quantity.  It  is  then 
also  necessary  to  introduce  a  Green  function  proportional  to 
«V;  £  and  factorization  now  must  be  assumed  for  quantities  of  the 
form  <ptpptpA>,  etc.  This  new  set  of  equations  propagates 
t ield-molecule  correlations  existing  in  the  initial  state,  but  does 
not  generate  new  correlations;  i.e.,  if  the  quantity  <P^PA>  factors 
at  the  initial  time,  it  factors  at  later  times  as  well.  Spontaneous 
emission  r f  radiation  is  not  generated  by  the  equations  for  this  set 
of  Green  functions,  and  one  must  go  to  the  next  higher  order  of 
accuracy  and  consider  a  Green  function  proportional  to  <P^PP^P>. 

The  set  of  equations  for  Green  functions  correct  to  this  order  is 
straightforward  to  derive,  but  too  cumbersome  to  illustrate  as  part  of 
this  publication. 

It  should  now  be  sufficiently  clear  how  the  technique  may  be 
applied  to  calculate  equations  describing  any  desired  physical  effect. 
Consideration  of  higher  order  effects  by  such  methods  usually  must  be 
carried  out  for  the  simplest  model  systems.  Nevertheless,  the 
hierarchy  of  Green  function  equations  does  allow  one  to  determine  in  a 
rather  clear  manner  exactly  how  various  physical  effects  are 
interrelated  and  what  the  effects  of  various  approximations  may  be. 
Even  with  neglect  of  the  intermolecular  interaction  term,  the 
semir lassical  Equations  (5)  still  contain  interesting  features.  If 
one  ignores  molecular  motion  by  dropping  the  differential  operator 
2  such  that  it  is  sufficient  to  consider  only  the  case  of 

X'  =  then  one  obtains  a  set  of  equations  equivalent  (for  a  finite 
set  of  molecular  levels)  to  the  Maxwell-Schrodinger  model.  Thus  even 
in  the  dipole  approximation  one  cannot  consider  a  molecule  as  a  point 
particle  in  the  semiclassical  approximation,  as  evident  by  the 
structure  of  the  kinetic  energy  term  in  Equations  (5),  which  cannot  be 


expressed  in  terms  of  a  single  center-of-mass  coordinate  and  internal 
molecular  variables. 

As  one  illustration  of  the  basic  theory  in  the  semiclassical 
approximation  (using  a  dipole  approximation  and  ignoring  the 
intermolecular  interaction  term),  we  have  considered  a  typical 
three-level  swept-gain  superradiance  model  in  which  relaxation  times 
and  diffraction  losses  are  introduced  empirically  and  an  initial 
tipping  angle  (Bloch  model)  replaces  the  Inclusion  of  spontaneous 
emission  [2,6].  In  this  example,  we  have  shown  that  the  growth  of 
cooperative  macroscopic  polarization  in  the  medium  during  the  pumping 
process  scales  as  the  ratio  of  the  pump's  temporal  duration  p  to 
the  characteristic  superradiance  time  r,  rather  than  to  D  (the 
delay  time),  as  has  been  speculated  previously  in  the  literature  [13]. 
The  influence  of  the  pump  pulse  shape,  as  well  as  area,  was  also 
determined  to  be  critical  with  regard  to  subsequent  cooperative 
relaxation.  Also,  two-photon,  coherent,  stimulated  Raman 
contributions  were  shown  to  make  important  contributions  to  pulse 
evolution  during  the  pump  time  regime.  This  finding,  in  turn,  has 
opened  up  a  new  area  of  research  and  investigation  which  we  are  now  in 
the  process  of  vigorously  pursuing.  Other,  significant  and  meaningful 
examples  to  illustrate  the  basic  theory  can  easily  be  constructed, 
while  important  effects  such  as  spontaneous  emission,  diffraction, 
etc.,  can  be  rigorously  treated  or  accounted  for  solely  upon  the 
inclusion  of  higher  order  levels  of  approximation. 


III.  AN  EFFECTIVE  HAMILTONIAN  FOR  THE  TREATMENT  OF  N 
IDENTICAL  TWO-LEVEL  ATOMS  INTERACTING  WITH  THE 
ELECTROMAGNETIC  FIELD 


Our  interests  in  the  Dicke  model  via  superradiance ,  optical 
bistability,  etc.,  have  led  us  to  examine  in  some  detail  the 
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interaction  of  a  set  of  N  Identical  atoms  with  the  electromagnetic 
field  [4,6].  In  effect,  the  field  creation  and  annihilation  operators 
conform  to  standard  Bose  commutation  relations  involving  two 
specific  modes,  k  and  £,  while  the  atomic  operators  Cj£*cik  are 
Bose  operators  [14]  which  satisfy  ti>e  commutator  relations 


(6) 


where  i,j  identifies  the  specific  atomic  level  and  k,£  identifies  the 
specific  atom.  In  this  way  one  attempts  to  treat  a  system  which  can 
be  identified  by  the  existence  of  a  many-particle  Hamiltonian  and 
single-particle  eigenstates,  thereby  assuring  the  existence  of 
elementary  anti-commutator  properties  for  the  products  C^C^p.  As 
a  consequence,  for  a  two-level  system  (energies  >  e^),  we  may 
choose  a  vector  operator  f^- , 


(C9  .C, 


2121 


p5*  C  +C 

c  L2ihi 


ci Tci i ) 


f l  *  2**  C.JC 


ir2i 


(7) 


and  transform  to  a  representation  obeying  the  commutation 
(-)/anticommutation  (+)  relations 


[fr  4  S^2si, 

m 


km 


(8) 


where  can  be  established  on  the  basis  of  Cj?>  and 

their  relations  in  a  system  involving  only  single  particle  eigen¬ 
states.  It  is  also  to  be  noted  that  the  three-dimensional  operator 

a  i  i  i 

f*  (components  f^,  f2»  f 3)  acts  as  an  angular  momentum 
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operator  in  which  fj  denotes  the  direction  of  diagonalization  or 
measurement.  If  H^-,  a  Hamiltonian  in  f_^»  has  the  form 

H1  *  r  •  f1  ,  Tj  s  E  cose  f  r2  =  2-,sEsine  e"1ct  , 

.  <  («> 
r3  *  Sine  e1a  , 

where  r_  is  not  explicitly  dependent  upon  time,  it  is  a  simple  matter 
to  rotate  f}  about  the  _6  axis  to  obtain  a  frame  of  reference  in 
which  H*  is  diagonal.  The  transformation  T  (expressed  in  terms  of 
Euler's  angles),  which  rotates  f^,  then  can  be  used  to 
evaluate  the  temporal  and  thermal  variation  of  fi  in  a  frame  of 
reference  in  which  the  generator  H*  is  not  diagonal,  assuming  simple 
relations  exist  for  this  variation  in  that  frame  of  reference  in  which 
H*  is  diagonal.  Cases  in  which  the  temporal  and  thermal  generators 
are  diagonal  in  different  frames  can  also  be  adequately  handled, 
requiring  (of  course)  the  explicit  use  of  more  than  one 
transformation,  when  three  different  frames  (i.e.,  fi(l),  fi(2), 
fi(3))  are  considered  in  the  mathematical  expression  of  one  physical 
quantity.  Our  initial  argument  is  then  that  if  one  can  reduce  a 
Hamiltonian  for  the  interaction  of  N  identical  two-level  atoms  with  an 
electromagnetic  field  to  the  form  specified  by  expression  (9),  then 
one  is  in  a  position  to  manipulate  it  with  relative  ease  and  rigor. 

Denoting  the  (internal)  field  by  ^{a^.a^land  the  atoms  by 
f_  ■  Z  f} ,  let  the  temporal  variation  of  a  general  operator  X  be 
established  in  this  environment  through  its  Hamiltonian  generator  Hx 
and  the  resulting  equation  of  motion.  For  complex  one  may  define 
also  its  transform 

X(w)  *  6x(w;Hx)  '  (10) 

as  an  analytic  two-branch  operator  defined  everywhere  outside  the  real 
axis  [15]  and  obeying  an  equation  of  motion 
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ttX(w)  *  (2ir)”^X  +  Z(w)  ,  Z 


(ll) 


As  our  system  spans  the  product  space  of  a_  and  £,  Hx  would  assume  a 
functional  representation 


and  the  temporal  variation  of  X  would  then  evolve  from  expressions 
(11)  and  (12)  and  the  statistics  for  a  and  An  inherent  problem  in 
this  configuration  is  the  hierarchy  of  equations  generated  from 
relation  (11)  when  one  handles  a  (or  the  electromagnetic  field)  in  its 
present  representation.  In  this  light  (as  a  first  step  toward 
achieving  expression  (9)),  it  would  be  desirable  to  project  £  into  the 
space  of  £  such  that  Equation  (11)  is  at  least  solvable  exactly  at  the 
single-particle  level  with  no  explicit  time  dependence  appearing  in 
the  Hamiltonian.  In  so  doing  we  must  abandon  the  statistics  of  £  and 
establish  a  mapping  of  £  onto  jF  (and  its  statistics)  that  is 
consistent  within  the  framework  of  the  theory  and  sufficiently  general 
to  warrant  its  application.  It  will,  of  course,  not  be  exact  or 
absolute  as  it  is  a  projection  or  approximation  that  in  practice  will 
depend  on  the  constraints  of  modeling  and  interpretation.  In  so  doing 
our  new  representation  of  £,  and  thus  the  total  system,  will  be  in  £ 
rather  than  £  ®  £,  and  it  should  not  be  looked  upon  as  an 
approximation  that  one  introduces  arbitrarily  or  at  random. 

In  order  to  effectively  stabilize  our  system  and  make  it 
accessible  to  statistical  interpretation,  we  shall  now  introduce  a 
phenomenological  damping  coefficient  (rate  of  photon  escape  from 
the  quantization  volume  V)  for  the  kc^  photon  mode.  This  device 
will  effectively  account  for  the  reflectivity  from  cavity  mirrors 
which  limit  the  quantization  volume,  etc.  and  will  emphasize  the 
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Markovian  nature  of  the  system  for  times  long  compared  with  some 
characteristic  relaxation  time  tc  *  The  electromagnetic 

field  we  (in  turn)  interpret  to  exhibit  both  internal  and  external 
origins,  with  the  latter  potentially  identifiable  as  an  applied  field 
in  a  coherent  state,  with  field  amplitude  a  and  carrier  frequency 
The  generator  Ha^,  we  shall  assume  to  induce  an  equation  of 
motion  for  a^  of  the  form 


i  3ak/3t  *  +  £  •  f  ,  y^  »  -  Imf^)  , 

with  a  similar  equation  for  the  conjugate  a£.  Taking  the 
of  Equation  (13),  using  expression  (11),  then  establishes 
requirement  that  if 

-►  0  ,  as  (w-fy)  -*•  0  ,  (ia) 

then  [16] 

\  *  -  2*  xk.  f(nk)a  , 

(15) 


(13) 

transform 

the 


with  a  similar  statement  for  aj^.  Thus  f_  acts  as  an  operator  whose 
temporal  variation  provides  a  means  for  transmitting  the  statistics  of 
£  (contained  in  Ha)  into  fCQ^)  and  its  conjugate,  and,  as  such. 
Equation  (15)  must  still  be  identified  as  existing  In  a  ®  f.  Require¬ 
ment  (1A),  on  the  other  hand,  implies  that  f^  is  capable  of  supporting 
£  to  the  extent  that  the  interaction  af_  will  prevent  ( fife )  from 

tending  to  become  singular  or  otherwise  unconstrained  in  its  integral 
over  time.  In  any  event,  the  field  £  has  now  been  correlated  to  the 
atoms  f_  through  the  interaction  (coefficients  Xj);  and  the  external 
field  (or  source)  may,  in  turn,  either  drive  the  atoms  or  be  coupled 
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r 


directly  to  the  Internal  field  £  through  the  coefficient  A 4  as  It 
appears  in  Equation  (IS). 

Substituting  Equation  (15)  into  expression  (12)  and  defining  [17] 

•  Ff  a  ^Fa  +  Fa+)  »  (16 

one  is  now  led  to  insert  this  expression  into  the  transform  (10) 
thereby  obtaining 


(17) 

In  general,  expression  (17),  as  an  operator  equation,  cannot  be 
directly  evaluated,  and  one  must  go  to  some  statistical  average  or 
measurement  in  order  to  extract  information  useful  for  physical 
interpretation.  The  statistical  average  of  the  product  of  ^fC'cfc)  and 
an  arbitrary  parameter  Q,  or  Q  >,  is  a  likely  basis  for 

evaluation.  One  is  now  in  a  position  to  reduce  expression  (17)  to  the 
space  of  f_  by  defining  the  respective  mappings 

<1(^)0  >  *  A(nk)  <fQ  >  ,  *  B(°|C)JF  »  (18) 


where  is  now  arbitrary  in  nature  and  _f  could  just  as  well  be 
represented  by  its  single-particle  version  f_i.  For  the  general  case 
of  many  particles,  relations  (18)  can  only  be  identified  in  the  realm 
of  approximations  to  the  system  and  the  mappings  (matrices)  A  and  B 
are  impacted  by  modeling  considerations.  For  restrictive  situations  A 
can  be  taken  identical  to  B  and  the  statistical  average  of  expression 
(17)  becomes  an  equation  for  self  consistency  in  B.  In  the  case  of  a 
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many-particle  system,  the  substitution  of  ftefc),  as  approximated  by 
relation  (18)  into  Hf,  as  defined  by  expression  (11),  does  not 
immediately  produce  H*,  as  defined  in  expression  (9).  One  must 
first  resort  to  a  "mean-field"  form  of  approximation  in  order  to 
decouple  operator  products  involving  different  atoms  (i.e.,  f*  f^, 
i  +  j).  Thus,  in  the  fabrication  of  a  single  atom  Hamiltonian  H*, 
the  operator  f^  must  be  replaced  by  its  statistical  average  <fj|>. 
The  mean-field  approximation  is  certainly  justified  in  models  which 
asstime  that  all  of  the  atoms  are  contained  in  a  volume  much  smaller 
than  an  atom  resonance  wavelength.  In  summary,  the  three  steps 
necessary  for  the  fabrication  of  an  H*  are:  (a)  the  correlation  of 
afc  to  fGlfe);  (b)  the  mapping  of  ffc^)  into  f^;  and  (c)  the 
factoring  of  f*  f£  in  a  mean-field  approximation.  Of  primary 
importance  is,  of  course,  the  choice  of  A  and  B  which  impacts  one's 
ability  to  successfully  model  any  systems  considered. 


Although  we  have  mentioned  the  existence  of  an  applied  external 
field,  we  have  as  yet  not  taken  into  account  the  impact  of  its 
potential  explicit  dependence  on  time.  Any  explicit  time  dependence 
must  be  removed  initially  prior  to  other  considerations,  otherwise 
fabrication  of  H*  (as  previously  envisioned)  is  jeopardized,  and 
statistical  and  equilibrium  properties  of  a  given  system  cannot  be 
effectively  treated.  A  convenient  solution  is  to  canonically 
transform  to  a  doubly  rotating  frame  of  reference  rotating  at  the 
carrier  frequency  and  discuss  statistics  and  thermodynamics  in 
the  rotating  frame.  In  this  frame  the  applied  electromagnetic  field 
appears  as  a  dc  field.  After  canonical  frequency  renormalization,  the 
transformed  Hamiltonian  Hj  will  assume  the  form  UHU”^  for  an 
appropriate  transformation  operator  U.  The  existence  of  a  given  U  may 
impose  certain  constraints  or  restrictions  on  the  generality  of  H, 
particularly  in  regard  to  the  preservation  of  counterrotating  terms 
which  may  defy  manipulation.  In  order  to  discuss  thermodynamic 
equilibrium  properties  in  the  presence  of  a  time-dependent  external 
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field,  one  must  introduce  an  effective  or  “spin"  temper* ture. 

Agreement  with  nonequilibrium  models  would  then  exist  in  the  limit  of 
zero  effective  temperature.  The  combination  of  effective  or  “spin" 
temperature  and  transformation  to  a  frame  rotating  at  the  carrier 
frequency  has  been  employed  effectively  by  two  of  us  (CMB  and 
CCS)  in  the  treatment  of  first-  and  second-order  phase  transitions  in 
the  extended  Dicke  model  [4]. 

IV.  LIGHT-ATOM  INTERACTION  IN  COOPERATIVE  PHENOMENA  IN 

QUANTUM  OPTICS 

As  introduced  in  the  previous  section,  one  of  the  principal 
theoretical  efforts  in  quantum  optics  in  our  group  is  concerned  with 
the  atom-atom  interaction  via  the  electromagnetic  field  -  an  area  that 
has  drawn  considerable  attention  in  the  literature  for  some  time.  A 
complete  review  on  this  subject  is  given  by  Milonni  and  Knight  [18], 

A  second  area  of  emphasis  and  study  in  our  program  is  an  investigation 
of  the  thermodynamic  properties  of  the  Dicke  model  [17-19],  defined  by 
the  Hamiltonian 


H  '^g(o1a  +  °fa+)  +3Peo1  •  (19) 

41  •  z  t 

where  a j  and  Oj  are  the  spin  operators  for  atom  i  and  a' (a) 

are  the  field  creation  (annihilation)  operators  for  the  resonant 

single  mode. 

It  was  first  shown  by  Hepp  and  Lieb  [19]  that  a  second-order  phase 
transition  appears  at  a  critical  temperature  Tc  if  huy'8Trd2p  <  1, 
where  d  is  the  atomic  dipole  moment  matrix  element,  w  the  frequency, 
and  p  *  N/V  the  density  of  atoms  in  a  cavity.  Despite  considerable 
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progress  made  by  Wang  and  Hioe  [20]  and  Gilmore  and  Bowden  [21]  In  the 
simplification  of  the  calculation  and  the  study  of  the  effect  of  the 
nonlinear  term  In  H,  it  was  still  not  clear  as  to  the  physical 
mechanism  which  causes  the  phase  transition.  Subsequently,  we  have 
used  various  different  formalisms  to  study  the  problem  and  it  is  now 
clear  from  our  work  [10,11]  that  the  phase  transition  of  the  atomic 
system  is  simply  caused  by  the  atom-atom  pair  interaction  induced  by 
the  radiation  field, 

Hmt  *  £ 5  °i  •  <20) 


where  g  %  g2/e. 

Equation  (20)  shows  that  the  atomic  interaction  in  the  Dicke  model 
is  an  infinitely  long  range  force,  since  g  has  no  spatial  dependence. 
This  conclusion  is  valid  for  a  single  mode  of  the  radiation  field. 
Otherwise,  we  must  assume  that  the  volume  Vg  containing  the  atoms  is 
much  smaller  than  the  resonant  wavelength,  which  therefore  is  much 
smaller  than  the  field  quantization  or  cavity  volume  V.  When  the 
inter-atomic  distance  becomes  larger  than  the  resonant  wavelength, 
then  g  becomes  equivalent  to  the  Van  der  Waal  force  between  two  atoms 
[18].  These  observations  seem  to  have  been  neglected  in  the 
literature.  We  also  have  extended  our  calculations  to  three-level 
atoms  in  the  Dicke  model  and  have  shown  the  existence  of  two 
second-order  phase  transitions. 

Optical  Bistability  (OB)  [22]  has  drawn  a  great  deal  of  attention 
in  the  last  few  years  after  the  pioneering  experimental  work  by  Gibbs, 
McCall,  and  Venkatesan,  and  McCall's  own  theoretical  effort  [23]. 
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Theoretical  models  [24,25]  for  the  analytical  interpretation  of  the 
hysteresis  of  the  input  field  Ex  versus  the  transmitted  field 
are  based  on  the  coupled  Maxwell-Bloch  equations.  The  theoretical 
results  obtained  show  that  the  relation  E-j  versus  Ex  is  quali¬ 
tatively  consistent  with  experimental  observations  for  the  absorptive 
case.  OB  based  on  this  model  [26]  is  essentially  a  saturation 
phenomenon  similar  to  magnetic  resonance,  and  (as  a  matter  of  fact) 
the  condition  for  OB  implies  specifically  that  the  atomic  absorption 
must  be  saturable. 

In  the  spirit  of  the  study  of  atom-field  interactions,  we  propose 
a  new  model  for  OB  which  is  assumed  to  be  represented  by  a  thermo¬ 
dynamic  ensemble  in  the  rotating  frame  of  the  external  field.  After 
transformation  to  the  rotating  frame  of  reference,  we  obtain  an 
atom-atom  interaction  similar  to  Equation  (20),  except  that 
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2  (mc-“o 
°  ("c 2  ♦  Y* 


(21) 


where  wc  is  the  cavity  frequency  of  the  Fabry-Perot  interferometer,  Y 
is  the  photon  escape  rate,  and 


a 


[(e  -  u0)2  +  A2] 


(22) 


is  the  renormalized  atomic  level  separation  with  A  to  be  deter¬ 
mined  self-consistently  from  the  equation  of  motion.  Distinctive 
characteristics  of  our  model  are  that  we  take  A  into  account  in  the 
operator  equation  for  a^(a)  and  study  OB  in  the  limit  of  large  Tj 
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and  T2  —  the  small  absorption  limit.  The  resultant  equation  for 
Ef  versus  Ej  is  then  given  by 


tanh  (o/2KT  ) 

h  ■  h  *  2Et5 - y— *- 


(23) 


where  Ts  is  the  effective  temperature  directly  analogous  to  the 
"spin  temperature"  in  the  rotating  frame,  and  A  (as  specified  in 
Equation  (22))  is  simply  Ej. 

In  the  limit  o)c  mw0  =  e,  the  condition  for  perfect  tuning, 
Equations  (21)  and  (22)  reduce  to  the  well-known  formula  (for 
establishing  OB)  given  in  the  previous  literature  [24,25]  with 
different  coupling  constant  C.  Although  our  work  may  be  used  to 
explain  the  dispersive  effect  better  than  the  previous  works  [24,25], 
we  should  emphasize  our  fundamental  difference  in  approach  and 
physical  mechanism  to  explain  OB.  It  seems  our  model  (by  using  the 
Green's  function  approach)  requires  and  T2  to  be  sufficiently 
large  so  that  each  atomic  level  is  not  critically  broadened;  whereas 
the  Maxwell-Bloch  equations  assume  small  and  T2  for  observation 
of  OB. 


Experimentally,  the  Stark  shift  A  (in  our  model)  can  easily  be 
observed  In  the  light  spectrum;  likewise  the  change  in  the  linewidth 
discussed  in  the  Maxwell-Bloch  model  can  also  be  observed  easily.  We 
hope  some  experimental  work  in  this  area  will  appear  soon  to  clarify 
the  situation.  Regardless  of  the  outcome  as  to  which  model  is  more 
adequate  to  describe  OB,  we  believe  it  interesting  that  a  first-order 
phase  transition  can  indeed  be  deduced  from  the  Dicke  model  when  an 
external  field  is  Introduced  to  modulate  the  atomic  interaction.  At 
present,  we  are  working  on  the  effect  of  introducing  finite  Tj  and 
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T2  into  our  model  and  including  the  spatial  dependetice  of  the 
interaction  discussed  earlier.  Likewise,  an  attempt  will  be  made  to 
extend  our  modeling  activities  to  larger  and  more  complex  systems. 


V .  CONCLUSION 


In  the  preceding  sections  we  have  outlined  the  scope  of  our 
program  in  cooperative  phenomena  and  matter/field  interactions  and 
detailed  some  specific  areas  of  current  thrust  and  momentum.  We  have 
tried  to  emphasize  the  development  of  appropriate  building  blocks  from 
scientific  first  principles  and  the  establishment  of  effective 
techniques  for  the  structuring  and  analysis  of  theory  and  data.  The 
theory  developed  in  Section  II  and  depicted  in  Equation  (5)  is 
exceedingly  general  and  definitive  both  in  its  scope  and  content.  Its 
application  to  the  case  of  swept-gain  superradiance  is  but  one 
illustration  of  its  general  appeal.  Even  here  a  rigorous  treatment  of 
effects,  such  as  relaxation,  diffraction,  and  spontaneous  emission,  is 
available  upon  consideration  of  higher  levels  of  detail  and 
approximation. 

Concurrently,  our  investigations  into  OB  have  characterized  that 
phenomenon  from  a  purely  thermodynamic  point  of  view.  As  a 
dipole-dipole  pair  correlation,  the  effective  interaction  is  identical 
to  the  retarded  dipole-dipole  interaction  that  gives  rise  to 
cooperative  radiation  reaction  and  frequency  shifts  in 
superfluorescence  [27]  and  is  the  same  effective  interaction 
responsible  for  the  existence  of  the  ’’superradiant"  phase  transition 
in  thermodynamic  equilibrium  [11,19].  Likewise,  the  Interaction  is 
also  identical  in  form  to  the  interaction  coupling  Cooper  pairs  in  the 
BCS  theory  of  superconductivity  [11,28],  Our  results  provide  not  only 
a  description  of  OB  as  a  first-order  phase  transition,  but  also,  give 
conditions  for  OB  in  terms  of  characteristic  material  and  cavity 
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parameters,  and  incorporate  the  effects  of  dispersion  as  well  as 
absorption.  In  Sections  III  and  IV  we  have  attempted  to  present 
certain  basic  aspects  of  this  research  that  impact  the  development  of 
the  formalism  and  provide  insight  into  its  structure  and  implications. 

In  future  research,  we  will  increase  the  rigor  and  complexity  of 
our  present  modeling  applications  and  generalize  our  theoretical 
activity  to  the  consideration  and  understanding  of  new  systems  and  new 
and  related  phenomena  and  effects.  Applications  involving  strong 
oscillating  fields,  multiple-level  systems,  a  combined  temporal  and 
spacial  variation,  the  rigorous  treatment  of  relaxation,  diffraction, 
spontaneous  emission,  etc.  will  be  considered  both  from  the  standpoint 
of  evaluation  and  implementation. 
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